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Abstract-In this paper, we study the unsteady 
hydromagnetic flow of a Walter's fluid (Model B') down 
an open inclined channel of width 2a and depth d 
under gravity, the walls of the channel being normal to 
the surface of the bottom under the influence of a 
uniform transverse magnetic field. A uniform 
tangential stress is applied at the free surface in the 
direction of flow. We have evaluated the velocity 
distribution by using Laplace transform and finite 
Fourier Sine transform technique. The velocity 
distribution has been obtained taking different form of 
time dependent pressure gradient g(t), viz., i) constant 
ii) exponential decreasing function of time and iii) 
Cosine function of time. The effects of magnetic 
parameter M, Reynolds number R and the viscoelastic 
parameter K are discussed on the velocity distribution 
in three different cases. 

Key words: Walter's B' fluid, open inclined channel, 
Laplace transform and finite Fourier Sine transform 
technique. 

1. INTRODUCTION 

A flowing liquid is said to have a free surface 
when the upper part of the bounding surface of the 
liquid is in contact with the overlying atmosphere, 
rather than with a solid, as would be the case of the 
flow were in a pipe completely full of the liquid. A 
flow with a free surface proceeding in a natural or 
artificial channel or conduit is called an open- 
channels flow. Open-channels may be divided into 
two types namely (i) Natural channels and (ii) 
Artificial channels. Natural channels range inform 
from the bounder-strewn bed of a mountain torrent to 
the relatively uniform channel of a large river. 
Artificial channels are man-made and are 
constructed in many forms. A pipe in which water 
flows with a free surface and flume of rectangular 



cross section constructed of sheet iron are two kinds 
of artificial channels. Other types are canals 
excavated in earth or blasted in rock, which either 
are left unlined or lined with smooth concrete or 
another suitable material. Unlined or lined tunnels 
bored through rock may also contain free surface 
flows. 

The flow of a liquid in an open inclined 
channel with a free surface has a wide application in 
the designs of drainage, irrigation canals, flood 
discharge channels and coating to paper rolls etc. 
Hence the flow of a liquid in an open inclined 
channel with a free surface under gravity has long 
been studied experimentally and several interesting 
empirical results have been reported by many 
investigators [3, 6, 7, 10, 11, 14]. The steady laminar 
flow of a viscous fluid flowing down an open inclined 
channel has been discussed by Satyaprakash [13], 
Gupta et al [4] have studied the flow of a viscous 
fluid through a porous medium down an open 
inclined channel. Venkataramana and Bathaiah [18] 
have studied the flow of a hydromagnetic viscous 
fluid down an open inclined channel with naturally 
permeable bed under the influence of a uniform 
transverse magnetic field. Unsteady laminar flow of 
an incompressible viscous fluid between porous, 
parallel flat plates has been investigated by Singh 
[12], taking (i) both plates are at rest and (ii) 
Generalized plane Coutte flow. The free surface was 
exposed to atmospheric pressure and bottom was 
taken as impermeable. Bakhmeteff [1], Henderson 
[5] and Chow [2] have discussed many types of open 
channel flows. Recently, many authors [8, 9 and 14] 
have studied the flow of Walter's B' fluid. 



The subject of Rheology is of great 
technological importance in many branches of 
industry. The problem arises of designing apparatus 
to transport or to process substances which cannot 
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be governed by the classical stress-strain velocity 
relations. Example of such substances and process 
are many, the extrusion of plastics, in the 
manufacture of rayon, nylon or other textile fibres, 
viscoelastic effects transported or forced through 
spinnerts and in the manufacture of lubricating 
grease and rubber. 

Non-Newtonian fluids have wide importance 
in the present day technology and industries; the 
Walter's fluid is one of such fluid. The model of 
Walter's B fluid is chosen for our study as it involves 
non-Newtonian parameter. The Cauchy stress tensor 
T in such a fluid is related to the motion in the 
following manner 

T = -PI+2?] 0 e-2K 0 j- (1.1) 

In this equation P is the pressure, / is the Identity 
tensor and the rate of strain tensor e is defined by 



2e = Vv+(Vv) 



(1.2) 



where v is the velocity vector, V is the gradient 

8 



operator and 



St 



denotes the convicted 



differentiation of a tensor quantity in relation to the 
material in motion. The convicted differentiation of 
the rate of strain tensor is given by 

— = — +vVe-e-Vv-(Vv) -e (1.3) 
St dt y } 

Here r/ 0 and K c . are, respectively, the limiting 

viscosity at small rate of shear and the short memory 
coefficient which are defined through 



7j 0 = [N{r)dT 



(1.4) 



and 



K 0 = [rN{T)dT 



(1.5) 



N(t) being the relaxation spectrum as introduced 

by Walter's [19, 20]. This idealized model is a valid 
approximation of Walter's fluid (model B) taking very 
short memory into account, so that terms involving 



[t"N(t)cIt, 



n>2 



(1.6) 



have been neglected. 



In addition to equation (1.1), the equation of motion 
and continuity are 



Ve=0 



p(v-Vv) = VT 



(1.7) 



(1.8) 



In this paper, we study the unsteady 
hydromagnetic flow of a Walter's fluid (model B) 
down an open inclined channel under gravity of width 
2a and depth d, the walls of the channel being 
normal to the surface of the bottom, under the 
influence of uniform transverse magnetic field. A 
uniform tangential stress is applied at the free 
surface in the direction of flow. We have evaluated 
the velocity distribution by using Laplace Transform 
and Finite Fourier Sine Transform techniques. Here 
it is assumed that (i) the fluid flows in the steady 

state for t<0, (ii) Unsteady state occurs at t>0 
and (iii) the unsteady motion is influenced by time 
dependent pressure gradient. The velocity 
distribution has been obtained in some particular 

cases i.e. when (i) g(t) = c* (ii) g(t) = c e~ bt and 

(\\\) g(t) = c cos bt, where b and c are constants. 

The effects of magnetic parameter M, Reynolds 
number R and viscoelastic parameter K are 
investigated on the velocity distribution in three 
different cases. 
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2. FORMULATION AND SOLUTION OF THE PROBLEM 



We consider the unsteady Hydromagnetic flow of a 
Walter's fluid (model B) down an open inclined 
channel of width 2a and depth d under gravity, the 
walls of the channel being normal to the surface of 
the bottom under the influence of uniform transverse 
magnetic field. A uniform tangential stress S is 
applied at the free surface. The bottom of the 

channel is taken at angle (3(Q< (3 <7t 1 2) with the 

horizontal. The x-axis is taken along central line in 
the direction of the flow at the free surface, y-axis 
along the depth of the channel and z-axis along 
width of the channel. A uniform magnetic field of 

^ = 0 

dx 



intensity H Q is introduced in y-direction. Therefore 

the velocity and the magnetic field are given by 

<7 = (u,0,0) and H = (0, H 0 ,0) . The fluid being 

slightly conducting, the magnetic Reynolds number 
is much less than unity, so that the induced magnetic 
field can be neglected in comparison with the applied 
magnetic field (Sparrow and Cess [15]). In the 
absence of any input electric field the equations of 
continuity and motion of the unsteady hydromagnetic 
Walter's fluid (model B) flowing down an open 
inclined channel at t>0 are given by 

(2.1) 



du dp . n (d 2 u d 2 u\ r _ | d 3 u d 3 u | 2rr2 

dt dX \dy dz J {dtdy dtdz ) 



0 = -^ + pgcos]3 (2.3) 

0 = -^ (2.4) 

dz 

Where p = density of the fluid 

g = acceleration due to gravity 
p = pressure 



jil = coefficient of viscosity 



<J = electrical conductivity of the fluid 



jU e = magnetic permeability 
K () = viscoelastic parameter 



The boundary conditions are 
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t<0; ll = u () 

t > 0; z = ±a,u= o 

n du 
7=0,//— =S 

y=d u=0 



(2.5) 



where u 0 is the initial velocity. 



we introduce the following non-dimensional quantities 

u - ul U, x = xl d, y - yl d, z - zl d, 
p=p/pU 2 , K* = K 0 / pd 2 , 5* = 5/ pU 2 



t = told 1 



(2.6) 



where the depth of the channel d is the characteristic length and the mean flow velocity U is the 
characteristic velocity. 

In view of the equation (2.6), the equations (2.1) to (2.3) reduce to (dropping the superscript *) 



du „dP R . 0 d 2 u d 2 u J ^ 
= -R- — +— sin /? + — — + — —-K 



dt 



dX F 



df dz L 



d 6 u d 6 u 
+ 



3.. A 



V 



dtd/ dtdz 2 



-Mu 



J 



(2.7) 



du 
dx 



= 0 



(2.8) 



Where M = OjU e z H 2 d 2 1 pv 



Magnetic parameter 



R=Udlv 



Reynolds number 



F=U 2 lgd 



Froude number 



K=K 0 /pd 



viscoelastic parameter 



The non-dimensional boundary conditions are 



f < 0; u = u 0 

t>0; z = l(=±al d),u 

y=0du/dy=SR 

y=lu=0 



= o 



(2.9) 
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3. METHOD OF SOLUTION 



Assuming 



-R— — +— sin /? 
dX F 



= g(t) at f>0 
= P at t<0 



(3-1) 



substituting z = (2lf/iT)-1 in equation (2.7) reduces to 



OU ,,01171011' 



ay 4/ df [dtdf 4/ dtdf 



.2 A 



d u 71 d u 
+ 



-Mu 



J 



(3-2) 



and the boundary conditions are reduced to 



t<0; u = u 0 
t>0\ f = 0,7r\ u 

y=0,^=SR 
7=1, u = 0 



= o 



(3.3) 



Now, since u 0 is the initial velocity i.e. at t < 0 , therefore taking g{t) — P in equation (3.2) 



2 00 



1 - COS fflt \ 



n 



P( cosh Cy \ 5i? sinhC(l-/) 
C 2 \ cosh C ) C cosh C 



sin/?/ 



(3.4) 



Where 



C 2 = 0 2 +M, 



TXT 

~27 



Now to solve equation (3.2) we take Laplace transform of equation (3.2) with respect to t (Sneddon [16]) 



u(y, f,s)= £ u(y, f,s)e st dt, s>0 



(3.5) 



we get 



d 2 u 7T 2 d 2 u (M + s)- 1 ( pKMcoshCy PKQ 2 SRMK sinh C(l-y) - 



+ 



dy M 2 df 2 l-Ks l-Ks 



u = 



V 



f + 

C z cosh C C 



+ 



C 



cosh C 



-#( 5 )-"b 



J 
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(3-6) 



Where 



g(s)=£g(t)e«dt 



On taking the finite Fourier sine transform the equation (3.6) with respect to f (Sneddon [16]) 



u*(y, N, s) = (u{y, f, s)sin Nfdf 



(3.7) 



We get 



2 * 



<L1 

3/ 



-H 2 u* = 



f 



V 



1 - cos nn 
N(l-Ks) 



J 



f 



V 



PKQ 2 -, , P P{l + MK)coshCy SR(l + MK)sinh Cy 



-g(s) + 

C C C cosh C C cosh C 



+ 



J 



Where H l = (7 + 

l-Ks 



Now, the boundary conditions are reduced to 



(3-8) 



J=0 



3u* 5i?(l-cos7V;r) 



3/ 

y=l, u* = 0 



SN 



(3-9) 



Integrating equation (3.8) under the boundary conditions (3.9), we get 



u* = 



^1-cos Nk^ 
{ N J 



PKQ 2 f cosh Hy ^ 8 



-1 



C 2 H 2 (\-Ks){ cosh H J H 2 (\-Ks) 



+ 



^cosh//j/ ^ 



V 



cosh // 



H 2 C 2 (\-Ks) 



^ cosh Hy ^ 



V 



cosh // 



V 



SC 



v 



Cos/?/// coshC/ SR sinh C '(1 -j) 
cosh // cosh C SC cosh C 



Now, inverting the finite Fourier sine transform as given by (Sneddon [16]) 



O oo 

/,s) = -£«*(y.Ms)sintf/ 

^ AA=1 



(3.10) 



In equation (3.10) we get 
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r v 2f,fl-cosA/p^ 



y 



// 2 C 2 (1-ife) 



Pift? 2 f cosh/// ^ g(s) ^ cosh/// 1 
C 2 H 2 (1- Ks) { cosh // J H 2 (l-Ks){ cosh // 

cos Hy A P ( cosh ify cosh Cj^i SR sinh C ( 1 - y) 



J 



v 



cosh// 



-1 



+ 



v 



cosh H cosh C 



J 



SC cosh C 
(3.11) 



sin Nf 



On inverting Laplace transform as defined by (Sneddon[16]) 



u{y,f,t) = —[ . u(y,f,s)e st dt 



In equation (3.1 1 ), we obtain 



2 I- cos NxX^, 4P(-l) r e A ' cos a r y r . 
*M N x(2r+l)(a r 2 + c 2 ) * ^ )8K ' 



5i? sinhC(l- 7 ) 1 . 
-J sin Nf 



C cosh C 



(3.12) 



Where 



r=0 



. . - 4(-l) r cosaye Ku 

h ( u ) = Y M T . a r =-(2r+l) 

V } ^7r(2r+l){l-K(a r 2 + Q 2 )} 2 V ' 



and 



4 = 



a r 2 + c 2 



\-k(a 2 + Q 



Particular cases 



Case 1. When g(t) = C* 



Using in equation (3.12), we get 



^_ 2yf l-cosAfrr ^4(-l) r 



ft N=l 



N 



r=0 



7c{2r+\)[ 



e ' cosa r y | C*^ coshCj ^ 



a r 2 + c 2 



V 



cosh C 



5i?sinhC(l-j) 
C cosh C 



] sin Nf 



(3.13) 



If we take the limit M — > 0, K — > 0, ^(7) = P= C* in equation (3.13) then we get the velocity distribution in 



the case of nonmagnetic and Newtonian fluid. In this case the velocity distribution is 
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2 yf l-cosTW r^ 



j 



C ( _ cosh Qy\ SR sinh Q(l - y) 



v 



coshQ ) Q coshQ 



sin Nf 



(3.14) 



This is in agreement with Sathyaprakash [13] 



Case II 



When g(t) = C*e bt , b> 0 , using in equation (3.12), we get 



u _ 2 j^ l-cosMr Vj^ 4P(-l) r cos a r ye Kt 



N=l\ iV J r 



=o ;r(2r+l)(a r 2 + C 2 ) 



- 4C*(-l) r cosa r ye- bt -e- At SR sinh C(\- y), . ATr 
+V — ^ —J sin Nf (3.15) 

to n ( 2r+ 1) [l - # ( a r 2 + (? 2 )] ( A r - b) C cosh C 



Case 



When ^(7) = C* cos bt , using in equation (3.12), we get 



u _ 2 j^ l-cosMr \j^ 4P(-l) r cos a r ye~ At 



xi^iK N ) 7t(2r+\)[a r 2 + C 



~\C*(-\) r cosa r y(A r cosbt+bsmbt- A r e At ) SR sinh C(l- y)i AT 

+ > — sin Nf 

h x(2r+l){l-K(a r 2 + Q 2 )}(A r 2 +b 2 ) C coshC 



(3.16) 



4. CONCLUSIONS 

Three distinct time dependent pressure gradient namely (i) g(t) = C* (ii) g(t) = C* e~^and (iii) 
g[t) = C* cos bt have been chosen to discuss the velocity profiles. Figures (1) to (3), (4) to (6), (7) to (9) are 
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drawn to investigate the effects of magnetic parameter M, time t and viscoelastic parameter K on the velocity 
distribution u in three different cases (i), (ii), (iii) respectively. In all the three cases we noticed that the velocity 
distribution increases with the increase in M or t where as it decreases with increase in K. Further we observe 
that in all the three cases the velocity profiles in nonmagnetic case are greater than with the magnetic case. 
Table (1) is drawn to bring out the effects of Reynolds number R on velocity distribution in three different 
cases. We noticed that in all the three cases the velocity distribution decreases with the increase in R. Further 
it is observed that in all three cases the velocity profiles in nonmagnetic case are less than with the magnetic 
case. 

GRAPHS 




y — > y » 

Fig. i u against y for different m Fig.2 u against y for different t 
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Table I 

Case I U against y for different R 



R 


M 


y=0 


0.2 


0.4 


0.6 


0.8 


1 


2 


1 


30.17 
57 


28.72 
465 


24.47 
127 


17.82 
51 


9.406 
355 


0 


2 


0 


14.51 
567 


13.83 
189 


11.80 

35 


8.622 
231 


4.568 
540 


0 


4 


1 


30.16 
235 


28.72 
465 


24.47 
127 


17.82 
51 


9.406 
355 


0 


4 


0 


14.50 

283 


13.83 
189 


11.80 

35 


80622 
231 


4.568 
54 


0 


6 


1 


30.14 
952 


28.72 
465 


24.47 
127 


17.82 
51 


9.406 
355 


0 


6 


0 


14.48 
990 


13.83 
189 


11.80 

35 


8.622 
231 


4.568 
54 


0 


8 


1 


30.13 
668 


28.72 
465 


24.47 
127 


17.82 
51 


9.406 
355 


0 


8 


0 


14.47 
716 


13.83 
189 


11.80 

35 


8.622 
231 


4.568 
54 


0 


1 
0 


1 


30.12 
385 


28.72 
465 


24.47 
127 


17.82 
51 


9.406 
35 


0 


1 
0 


0 


14.46 
432 


13.83 
189 


11.80 

35 


8.622 
231 


4.568 
54 


0 



Table II 

Case II U against y for different R 



R 


M 


y=0 


0.2 


0.4 


0.6 


0.8 


1 


2 


1 


34.14 
086 


32.50 
528 


27.71 
082 


20.20 
846 


10.68 
239 


0 


2 


0 


17.33 
922 


16.52 
671 


14.12 

038 


10.33 
622 


5.493 
201 


0 


4 


1 


34.12 
803 


32.50 
528 


27.71 
082 


20.20 
846 


10.68 
239 


0 


4 


0 


17.32 
638 


16.52 
671 


14.12 

038 


10.33 
622 


5.493 
201 


0 


6 


1 


34.11 
519 


32.50 
528 


27.71 
082 


20.20 
846 


10.68 
239 


0 


6 


0 


17.31 
35 


16.52 
671 


14.12 

038 


10.33 
622 


5.493 
201 


0 


8 


1 


34.10 
235 


32.50 
528 


27.71 
082 


20.20 
846 


10.68 
239 


0 


8 


0 


17.30 
071 


16.52 
671 


14.12 

038 


10.33 
622 


5.493 
201 


0 


1 
0 


1 


34.08 
952 


32.50 
528 


27.71 
082 


20.20 
846 


10.68 
239 


0 


1 
0 


0 


17.28 
2788 


16.52 
671 


14.12 

038 


10.33 
622 


5.493 
201 


0 
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Table 



Case 



U against y for different R 



R 


M 


y=0 


0.2 


0.4 


0.6 


0.8 


1 


2 


1 


31.2881 


29.79183 


25.40008 


18.52646 


9.795776 


0 


2 


0 


15.78622 


15.04873 


12.86061 


9.417838 


5.008056 


0 


4 


1 


31.27597 


29.79183 


25.40008 


18.52646 


9.795776 


0 


4 


0 


15.77338 


15.04873 


12.86061 


9.417838 


5.008056 


0 


6 


1 


31.26314 


29.79183 


25.40008 


18.52646 


9.795776 


0 


6 


0 


15.76055 


15.04873 


12.86061 


9.417838 


5.008056 


0 


8 


1 


31.2503 


29.79183 


25.40008 


18.52646 


9.795776 


0 


8 


0 


15.74771 


15.04873 


12.86061 


9.417838 


5.008056 


0 


10 


1 


31.23747 


29.79183 


25.40008 


18.52646 


9.795776 


0 


10 


0 


15.73487 


15.04873 


12.86061 


9.417838 


5.008056 


0 
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